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Abstract 



Context. Outbursting AM CVn stars exhibit outbursts similar to those observed in different types of dwarf novae. Their light-curves 
combine the characteristic features of SU UMa, ER UMa, Z Cam, and WZ Sge-type systems but also show a variety of properties 
never observed in dwarf novae. The compactness of AM CVn orbits and their unusual chemical composition make these systems 
valuable testbeds for outburst models. 

Aims. We aim for a better understanding of the role of helium in the accretion disc instability mechanism, testing the model for dwarf 
novae outbursts in the case of AM CVn stars, and aim to explain the outburst light-curves of these ultra-compact binaries. 
Methods. We calculated the properties of the hydrogen-free AM CVn stars using our previously developed numerical code adapted 
to the different chemical composition of these systems and supplemented with formulae accounting for mass transfer rate variations, 
additional sources of the disc heating, and the primary's magnetic field. 

Results. We discovered how helium-dominated discs react to the thermal-viscous instability and were able to reproduce various 
features of the outburst cycles in the light-curves of AM CVn stars. 

Conclusions. The AM CVn outbursts can be explained by the suitably adapted dwarf-nova disc instability model but, as in the 
case of its application to hydrogen-dominated cataclysmic variables, one has to resort to additional mechanisms to account for the 
observed superoutbursts, dips, cycling states, and standstills. We show that the enhanced mass-transfer rate, due presumably to variable 
irradiation of the secondary, must not only be taken into account but is a determining factor that shapes AM CVn star outbursts. The 
cause of the variable secondary's irradiation has yet to be understood; the best candidate is the precession of a tilted/warped disc. 

Key words. Accretion, accretion discs - Instabilities - Stars: dwarf novae - binaries: close 



1. Introduction 

AM CVn stars (AM CVns) are binary systems with very short 
orbital period (10 - 65 min) in which a He or C/O white 
dwarf primary accretes matter lost by a Roche-lobe filling sec- 
o ndary star, whic h is believed to be also a helium white dwarf 
R Nelema nsl2005h . The observed optical spectra show that accre- 
tion discs in these systems are helium-dominated and hydrogen- 
free. This spectral feature is indeed the criterion for including an 
ultra-compact binary to the AM CVn class. AM CVns are ob- 
serve d in three distinct luminosity states (see e.g. lBildsten et al.l 
2006). Two of these are the persistent-luminous (high) and the 
outbursting states, which are similar to those seen in cataclysmic 
variable stars (CVs) in which the secondary is hydrogen-rich and 
usually non-degenerate. The third state, during which the system 
is persistent and faint, does not have an equivalent among CVs. 
Smak i(|1983l) remarked that the disc instability model (DIM; see 
Lasotal200lL for a review and references) which describes dwarf 



novae (DN) also applies to AM CVn stars and can explain the 
properties of AM CVn itself (high state) and of GP Com (faint 
state). Because the instability in question is triggered by opac- 
ity changes, the differences between AM CVns and CVs are 
supposed to be explained mainly by the higher ioni zation po- 
tential of helium compared to that of hydrogen (see ICannizzol 
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1 The nature of the two shortest orbital period, (X-ray pulsating) sys- 
tems HM Cnc and V407 Vul remains a mystery. 



11984 lTsugawa~& Osaki 1997; Las ota et alJ l2008). However, as 
we show in the present paper, other effects related to the helium 
atomic structure as well as the presence of metals may also in- 
fluence the outburst cycle properties. 

In this context it is also worth stressing that angular momen- 
tum transport in discs depends on the microphysics and may 
be significantly affected by the chemical composition (e.g. the 
chemical composition influences the Prandtl number, which af- 
fects the development of the magneto-rotational instability). One 
of the methods to test the disc angular momentum transport 
mechanisms is outburst modelling. 

Although still poor compared to that of CVs, the AM CVns 
observation database has been enriched by recent observations 
and we know more about thes e binaries than we did in 1997. 
Then, in their pioneering work.lTsugawa & Osaki (119971) used a 
highly simplified description of the DI M (Osaki 1989) and not 
the full DIM model (such as e.g. llchikawa & Osaki 19921) . 

These facts as well as the desire to ap ply our version of 
the DIM for hel ium-rich accretion discs (lLasota et al.l [2008; 
iKotko et alfeOlOl) to AM CVn outbursts has been the main mo- 
tivation for writing the present article. In Section |2] we discuss 
the general properties of AM CVn outbursts, concentrating on 
those we are trying to explain in the following parts of the article. 
Section|3]describes the DIM in the context of helium-dominated 
accretion discs. Section|4]briefly recalls what is known about the 
chemical composition of AM CVn discs and how they comply 
with the stability criteria set by the DIM. In Section|5]we inves- 
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tigate the properties of helium-rich disc models in general and in 
Sections[6]and[7]we apply them to various types of AM CVn star 
outbursts. We discuss our results and summarize the conclusions 



2. Outburst properties of AM CVn stars 

From among 27 systems classified as A M CVn binaries, 11 
have been observed to produce outbursts dLevitan et al.ll20lTat 
Ramsav e t alJl2012h . Because of their relative faintness (~ 14 - 
20 mag) the light-curves of these systems are not as well covered 
and sampled as those of the longer period, hydrogen-dominated 
CVs. Based on a very incomplete (but soon to be increased by 
observational camp aigns such as the Palomar Transient Factory, 
iLevitan et alj d201 lbl) ) AM CVn database one can conclude that 
most (all ?) outbursting helium-rich systems show superout- 
bursts or, to be more precise, outbursts that resemble superout- 
bursts observed in h ydrogen-rich dwarf novae such as SU UMa 
stars (Warner 2003). The most important similarity is the pres- 
ence of superhumps - photometric humps with periods slightly 
different from that of the orbital modulation. The amplitudes of 
these outbursts range from 3.5 to 6 mag; the recurrence times 
from ~ 45 to ~ 450 days, but this latter parameter has been deter- 
mined only for five systems with orbital periods between 24.52 
and 28.32 min. AM CVn with longer orbital periods are believed 
(e.g. ILevitan et al.ll201 lal) to have recurrence times longer than 
three months but this still requires confirmation. One should 
also stress that in two cases a substantial variation of the re- 
curre nce time has been observed: in CR Boo from 46.3 to 14.7 
days dKato et alj|2001h and in KL Dra from ~ 65 to ~ 44 days 
dRamsav et alJ I2010L 12012b . The superoutbursts of AM CVn 
stars typically last about 20 days but can b e as short as 9 days . 
On th e other hand, in some cases (see e.g. lPatterson et al]|1997l 
2000) superoutbursts are prolonged by a series of low- amplitude, 
frequent outbursts: the so-called cycling state that may last as 
long as the "smooth" superoutburst itself. Often the superout- 
burst is in fa ct not so smooth because it contains a fairly pro- 
nounced dip dRamsav et alJl2010U 2012). The relation, if any, of 
these dips to the cycling state and of these two features to the 
"dipping state" observed in the light-curves of some of WZ Sge- 
type stars has not been satisfactorily elucidated yet. 

Much less in known about the so-called "normal outbursts" 
that are common in hydrogen-rich dwarf novae, but rarely ob- 
served in outbursting AM CVn stars. 

From the observational point of view, "normal" outbursts are 
not very well defined but they are supposed to correspond to the 
quite narrow eruptions observed in U Gem-type dwarf novae. 
They are defined in opposition to the superoutbursts whose most 
characteristic feature are superhumps. Normal outbursts could 
therefore be defined by the lack of superhumps (but in some 
cases normal outbursts following a superoutburst do show su- 
perhumps.) 

In the case of CVs, no short period (with orbital periods be- 
low the 2 - 3 hr "period gap") dwarf nova shows a pure normal 
outburst cycle. The SU UMa stars show a supercycle in which 
normal outbursts appear between superoutbursts, and binaries of 
the WZ Sge-type show superoutbursts only. The shortness of pe- 
riod seems to be the decisive factor determining the presence of 
superoutbursts^ which is confirmed by the fact that all AM CVns 
exhibit these types of outbursts. 



Noticing that in CR Boo and V803 Cen the outbursts seen 
during the cycling state satisfy the empirical Kukarkin-Parenago 
relation, connecting the ampli tude A„ with the r ecurrence time 
T„ (see iKotko & Lasotl I2012L and Sect. [67Tb . Patt erson etal] 
< fl997ll2OO0F concluded that t hey can be qu alified as normal out- 
bursts. On the other hand. lKato et alj d2001l) considered this state 
as equivalent to the standstill in Z Cam type dwarf novae. Until 
now the only clear case for relatively well sampled normal out- 
bursts has been provi ded by observations of PTF1 J0719+4858 
dLevitan et alJ 1201 "Tab . In this case the amplitudes are ~ 2.5 
mag, the duration ~ 1 day and the recurrence time about 10 
days. In this system normal outbursts appear bet ween superout- 
bursts as in SU UMa stars. ILevitan et all (1201 lal) suggested that 
AM CVn with orbital periods shorter than ~ 27 min are equiv- 
alent to SU UMa stars while helium-rich systems with longer 
periods would correspond to long recurrence time WZ Sge stars 
showing only superoutbursts. Although not implausible, this hy- 
pothesis has to await a more complete observational coverage 
of AM CVn binaries. Let us just mention that a complete anal- 
ogy between dwarf novae and their helium-rich cousins might be 
difficult to establish since e.g. the short-period KL Dra seems to 
exhib it only (short-recurr ence period) superoutbursts (see, how- 
ever, iRamsav et alj|2012h . 



2.1. The DM for AM CVn outbursts: preliminary 
considerations 

A model accounting for the properties of the three categories 
of AM CVn stars - permanently bright, permanently faint, and 
outbursting has to take into account the two main differences 
between these systems and the hydrogen-dominated CVs: the 
helium-dominated accretion discs and shorter orbital periods, 
that is, short er (less extended) accretion discs. As noticed by 
ISmakl d 1 983b the obvious candidate here is the DIM that is used 
to explain outbursts of hydrogen-rich dwarf novae. It might seem 
that the adaptation of this model to the case of AM CVn is pretty 
straightforward because it should consist only in rescaling to 
helium-rich discs the critical values defining the instability strip. 
This simplistic view assumes, however, that the DIM is fully suc- 
cessful in explaining most of the dwarf nova outburst properties. 
Unfortunately this is not the case. 

First, the standard DIM does not account for superout- 
bursts, which is particularly embarrassing for AM CVn stars in 
which this type of outburst is observed most frequently. Two 
modifications of the DIM have been propose d to acco unt for 
super outbursts of dwarf novae. Osakil (1 1989b (see also |Osaki 
1996, for a review) proposed that superoutbursts in SU UMa 
stars are produced by a tidal-thermal instability (TTI) trig- 
gered by a large enhancement of the tidal torque acting on 
the disc. Despite successes in describing some of the superout- 
burst properties, the TTI model has encountered some insur- 
mountabl e difficulties when trying to ag ree with observations 



2 How short is the required period is subject to de bate since U Gem , 
at about 4 hr orbital period, showed a superoutburst dMason et al J 19881 : 
ISmak & Wa aeen 2004; Smak 2006; Schreiber 20071) 



mountabie dimcu mes wnen trying to ag ree witn observations 
(see e.g. iBuat-Menard & Hameurvll2002HSchreiber et al.ll2004t 
Hameurv & Lasotall2005l and referen c es ther ein). Recently, in 
a remarkable series of articles Smak (2OO9a,b]00) has under- 
mined the basis of the TTI model by casting doubt on the 
presence of an eccentric disc in the SU UMa stars and pro- 
viding strong observational arguments in favour of a different 
model: th e enhanced ma ss-transfer (EMT) m odel, originally pro - 
posed bylOsakil (11985b and deve l oped by lLasota et alj (U 995). 
Buat-Menard & Hameurvl d2002b . ISchreiber et al.l (12004 "and 
Hameurv & Lasotal (l2005h . According to the EMT model: su- 
peroutbursts are due and begin with a major enhancement in the 
mass transfer rate. During the "Hat-top" part of the superout- 
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burst the mass transfer rate decreases slowly, causing the ob- 
served luminosity to decline. The superoutburst ends when the 
mass transfer rate decreases below its critical value, result ing in 
a tran sition to the quiescent state of the dwarf nova cycle (Smak 
2008). In this article we will be testing the EMT superoutburst 
model for several reasons. 

First, there is (mostly indirect) evidence of large variations of 
the mass-t ransfer rate at l east in some outburstin g AM CVn (see 
e.g. Patterson et al. 2000). Second, as shown by ISchreiber et alj 
(2004), the differences between standard SU UMa stars and 
ER UMa systems are a natural outcome of the EMT, but not of 
the TTI model therefore it is unlikely that the latter will be able to 
describe a system behaving like a WZ Sge, S U UMa, ER UMa, 
and Z Cam star. Third, since (unchallenged) |ma3 (|2009a) has 
undermined the very basis of the TTI model for the SU UMa 
stars it is normal to assume it is also not relevant to AM CVn 
star outbursts. 

Although normal outbursts seem to be infrequent in 
AM CVn stars, they are the phenomena the DIM is supposed 
to describe and explain. Since the standard DIM can be tested 
only on such outbursts, there is an obvious difficulty in applying 
this model to AM CVns. The recent observation of AM CVn star 
outbursts that can be clearly identified as normal (Levita n et alj 
1201 lal) promises a more straightforward comparison with the 
DIM predictions. In particular, it might help solving one im- 
portant problem with this model: the scali ng and re-sca ling of 
the viscosity parameter a. As first noted by Smak] d 19841) . to ob- 
tain dwarf nova outbursts with observed amplitude the viscosity 
parameter during outburst must be at least four times larger 
than the quiescent ff c . When implemented, this ad hoc assump- 
tion increases the ratio of critical surface densities which un turn 
allows (because of the S-shaped equilibrium curves) adequately 
large outburst amp litudes . 

As showed by Smak (1999) for dwarf novae, ah ~ 0.2. He 
arrived at this conclusion in two ways. First, by comparing the 
observed relation between the normal-outburst decay time and 
the orbital period (fdec - ^orb) with the model-calculated decay 
time fdec = Rd/vr, where Rp is the outer disc radius and vr 
the viscous speed in the hot disc. Second, by using the normal- 
outburst width vs orbital period (W - P OI b) relation to compare 
with the width of the normal outbursts obtained from the model. 
The viscous speed is 



(1) 



use the two fiducial values - 0.1 and = 0.2, hoping that 
observations will soon bring more information. The ratio ah/a c 
required to reproduce observed properties of AM CVn outbursts 
have to be determined through the DIM calculations. This prob- 
lem will be addressed below. First, we have to establish how the 
change of chemical composition affects the DIM. 

3. The disc instability model for AM CVn stars: 
application 

The model assumes that the disc is geometrically thin, which al- 
lows decoupling the radial and vertical structure equations. The 
radial structure of the disc is described by the equations of mass, 
angul ar momentum and en ergy conservation (see Eqs. (1), (2), 
(4) in Hameur y et al.lll998T hereafter H98) and its vertical struc- 
ture is assumed to be in hydrostatic equilibrium, which is de- 
scribed by equations (14-16) and (23) from H98. 

By solving the equations of the disc local vertical structure 
one obtains the values of effective temperature T e ff (or central 
temperature T c ) and surface density E, for which the thermal 
equilibrium, at a given disc radius R, is calculated. At each R 
these solutions form on the T - £ plane the well-known S-shaped 
curve. The lower branch of this curve, terminating at a maxi- 
mum value of the surface density S~ ri , represents stable cold disc 
equilibria, while the upper branch starting at a critical minimum 
value of the surface density £+ rit , corresponds to stable, hot equi- 
librium disc solutions. Solutions on the middle branch are ther- 
mally unstable. Their existence is the cause of the dwarf nova 
outbursts. 

For a pure helium disc (7=1) the critical surface densities 
(S^ rit ) and the critical effective temperatures (T* B ) can be fitted 
by the formulae found by Lasota et al. (2008): 



eff 



528 <} 81 ^| 07 Mf 36 gcm- 2 



M0 



13000 < 01 C< 3 K 



S;„ = 1620 a-™*R]^M; 0M gem" 2 



r eff = 9700a o r*7o 



0.1 "10 J "i 
-0.01 d-0.09 ^0.03 



K, 



(2) 
(3) 
(4) 
(5) 



v 2 _j ykT c 1/2 

vr ~ - ~ a h c s v K ~ a h R 1 , 

R prni H 

where v is the kinematic viscosity coefficient, c s the sound speed, 
vk the Keplerian velocity, // the mean molecular weight, y the 
adiabatic index, m# the hydrogen atomic mass, and T c is the 
mid-plane temperature. Therefore the viscous speed in AM CVn 
discs is half that of hydrogen-dominated CVs (the critical tem- 
peratures, for which ionization/recombination of the dominant 
element in the disc takes place, are twice as high but the molec- 
ular weight is four times higher). Since AM CVn discs are 
shorter (by a factor ~ 2) than those of CVs, decay times in 
the former should be comparable to those in the latter. The 
available limited sampl e seems to confirm that. The correlation 
found by ISmakl d 1999b between the observed outburst widths 
W(d) (in days) and the orbital period P or b (in hours) is W(d) = 
(2.01 + 0.29)P orb (hr) (a78±011) . Therefore for P orb = 25 min the 
normal outburst duratio n should be ~ 1, day as observed in 
PTF 1J0719+4858 bv lLevitan etall (l2011al) . Based on this very 
limited evidence, one can conclude that the value of the viscosity 
parameter in h ot AM CV discs is clo se to that of hydrogen-rich 
dwarf novae dKofko&Lasotall2012h . For convenience we will vS ~ crT, 



whereas for a helium disc with solar metal abundance (Z = 0.02, 
Y = 0.98) one has 



K nt = 380 < 78 <? 6 Mr 035 gem- 2 
r e + ff = llSOOaor'^o - 1 



-0.08 M 0.03 K 



V- *n -0.82 nl lO j^-0.37 „ -2 

S crit = 612 a 0.l R 10 M l 8 cm 

T eff = 8690ao r< 09 <° 3 K, 



(6) 
(7) 
(8) 
(9) 



where ao.i is the viscosity parameter in units of 0.1, R\q is the 
disc radius in 10 10 cm and M\ the white dwarf primary's mass 
in solar units (see Appendix lAl for the complete list of critical 
values, also for Z = 0.04). 

The critical effective temperatures are independent of the vis- 
cosity parameter a (the very weak dependence in some cases re- 
sults from the imprecision of the fit) because in thermal equilib- 
rium the flux (h ence T e ff) is indep endent of the viscosity mech- 
anism (see e.g. iFrank et al.ll2~002h . This allows a fairly simple 
understanding of the S -curve properties. 

From mass and angular-momentum conservation equations 
one obtains 



eff- 



(10) 
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The kinematic viscosity coefficient is v = 2/3 a c s H, where c s is 
the speed of sound and H the disc semi-scale-height. Therefore 
(for a gas-pressure-dominated disc) one has 



aT r Y. ~ crT, 



elf' 



(ID 



which implies that for T c oc a" and E oc a m , n + m = — 1, Eq. (fTTb 
and the energy transfer equation linking T c and T^, determine 
the slopes of the S-curve's branches T e s{Y). 

In the case of radiative cooling the effective and central tem- 
peratures are related through 



T, = 



' 3 M 1/4 r 



(12) 



where T tot is the total disc optical depth (see iDubus et alj 
[1991 . For a fully ionized disc the Rosseland opacity coeffi- 
cient can be written as k ~ ~£H~ 1 T~~ 1 ^ 2 cm 2 /g from which fol- 
lows the well-known relation T e s oc £ 5 / 14 found in the classical 
Shakura & Sunvaevl (1 19731) solution. On the other hand, a very 
steep decrease of opacity with temperature induces a change of 
slope of the T e s(X) relation. This is what happens during recom- 
bination. 

Fig. Q] shows S-curves of discs with four chemical compo- 
sitions: 7 = 1.0,7 = 0.98,7 = 0.96, and solar abundance 
(X = 0.7, 7 = 0.28, Z = 0.02). As expected, for all cases the 
upper branches are parallel (and very close). Since the change 
of slope at Elj is caused by a recombination-induced change in 
the opacity dependence on temperature, it follows that the upper 
S-curve bend will be the highest for the element with the highest 
ionization potential. Hence £+ rft (7 = 1) > 2+ rft (7 = 0.98) > 
S^ rit (7 = 0.96) > Ei- t (solar). Cooling by convection, which 
flattens the temperature profile, affects the position of the up- 
per bend of the S-curve, but not the relative order of the critical 
surface-densities values. 

For the lower bend of the S -curve, the situation is sl ightly 
more complicated. According to Tsuga wa & Osakil (1 19971) it ap- 
pears where the disc becomes optically thin (r ~ 1) and they 
explain the higher value of E".. for helium by its low opacity 
that must accordingly be compensated for by higher surface den- 
sity to give kL ~ 1 . However, as seen from Table [1] this expla- 
nation cannot be correct. A s already noticed and explained by 
ICannizzo~& Wheeled ( fl98i see also ISmakl (119991) ) the optical 
depth at E~ rit is always (for the standard assumption a < 1) larger 
than 1 and the turnover at E~ ri . is related to the optical depth in a 
less straightforward way. 

From Table Q] one can see that at critical points the rela- 
tion Eq. ( fT2l is approximately satisfied only for 7=1. For 
the other two cases the temperature gradient is much flatter 
than predicted by Eq. (fT2l . reflecting the importance of con- 
vective energy transport, consistent with the high optical depths 
of these configurations. Therefore for 7=1 (low opacity) the 
turnover at S~ ri results from the change in opacity temperature- 
dependence, whereas for the high optical depth cases of Z = 0.02 
and Z = 0.04 the change of slope results from the (opacity- 
related) change in importance of convective energy transport. 

Both (maximum and minimum) critical surface-density val- 
ues are higher in helium dominated discs than in discs where 
hydro gen dominates. However, as has been realized by Smak 
(1984J), the S-curves used in the DIM must be modified if the 
model is to reproduce the outburst properties of dwarf novae. 
The viscosity parameter a cannot be constant during outburst, it 
must be larger in hot discs than in cold, quiescent discs. The DIM 
uses "effective" S-curves by combining configurations with two 
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Figure 1. S-curves for M\ = 1.0, R = 10' cm , a = 0.1, and four 
chemical compositions. 
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Figure 2. "Effective" S-curves for combined cold and hot branches 
corresponding to two values of a: a c = 0.05 and a h = 0.2. The curve 
on the left corresponds to solar composition, while that on the right 
corresponds to Y = 0.98. M, = 1.0, R = 1.0 X 10 9 cm. 
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108.0 


50.2 


39.0 
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10700 


11700 


11200 


6960 


W] 


15600 


14250 


14200 


8210 


K(Z,77)[cm 2 /g] 


4.2 x 10-' 


4.5 x 10 2 


5.4 x 10 2 


1.3 x 10 3 


Ttot(S, 77) 


45.3 


2.24 x 10 4 


2.11 x 10 4 


7.03 x 10 3 



Table 1. S-curve critical points for four abundances. R = 1.0 x 10 9 cm, 

a = 0.1 



different values of a. Observations of dwarf nova normal out- 
bursts imply that the ratio between the two is between 4 and 10. 
Examples of effective S-curves are shown in Fig. |2]for an a ra- 
tio equal 4.0. This ratio has been assumed to be the same for 
a hydrogen- and helium-dominated disc. Since the a variation 
used is purely arbitrary, gauged through observed dwarf nova 
light curves, this assumption is not necessarily valid. We will 
discuss this problem in more detail in Sect. 
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4. Mass transfer rates and stability 

The first to consider t he stability of helium accretion discs in 
AM CVn systems was ICannizzol (1984). Here we use the latest 
Opal opacity tables and a larger statistics of AM CVn stars to 
test the stability criteria versus the observed properties of these 
systems. As discussed above, the critical parameters defining the 
disc's stability depend on the chemical composition. 

4. 1 . Chemical composition of accretion discs in AM CVn 
stars. 

In general, observations clearly show that discs in AM CVn stars 
are not made of pure helium. We briefly describe the observed 
AM CVn spectra to show what their disc metallicity might be. 

The two persistently faint systems - GP Com and V396 Hya 
- are very special. Unlike the other (outbursting or persistently 
bright) AM CVn stars, they show a significant overabundance of 
NV, in addition to a strong Hel and a weaker Hell line, but not 
the Si line, that is usually seen in other binaries of this type. This 
suggests metal poor secondaries. Other binaries in a low state, 
such as SDSS J0902 or SDSS J1552 also show weak Sill and 
Fell lines. 

In the permanent high-state systems AM CVn or HP Lib the 
spectrum is dominated by the absorption lines, but sometimes 
emission in Hell is also detected. In the optical part of the spec- 
trum the asymmetrical, broad absorption line of Hel dominates. 
The UV spectru m shows broad a bsorption lines of Hell, NV, 
NIV, SilV, CIV (IWade et alJl2007l) . The spectra of the erupting 
systems, CR Boo and V803 Cen have similar characteristics dur- 
ing outbursts. 

For GP Com IStrohmaverl (I2004 estimated X He = 0.99 
and Z = 0.01. Because GP Com and V396 Hya are the two 
metal -poor systems, we can consider Z = 0.01 as lower limit 
for an AM CVn disc metallicity. However, the exact metallic - 
ities of particula r systems have not been well determined yet 
dNelemans et al.ll201Tl and private communication). 

Throughout, our fiducial model of outbursting AM CVn stars 
assumes Y = 0.98, Z = 0.02, but when discussing the general 
properties of the systems and the model we will also consider 
other possibilities. 

4.2. Stability of observed systems 

In the standard version of the DIM the mass transfer rate M a 
from the secondary is assumed to be constant. Its value deter- 
mines the stability of the disc. 

For pure helium dis cs the critical values of the accretion rate 
are (lLasota et al.ll2008l) 



MX = 1.01 X 10 17 ao°r 05 ^ 68M r a89 g^' 



and 



M~ = 3.17 x 10% ( [ 02 /4 66 Mr (X89 gs" 1 - 



(13) 



(14) 



To be in a hot (cold) stable equilibrium the accretion rate in a 
disc must be higher (lower) everywhere than the corresponding 
critical M. Therefore the stability conditions are 

- a stationary (M(R) = const.) disc is hot and stable when 

M tr >M+(R D );, 

- a stationary (M(R) = const.) disc is cold and stable when 
M tt <M a (R ln ). 
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Figure 3. The log M u - P mb plane. Dots - persistent AM CVn systems 
with known distances; asterisks - outbursting AM CVn systems with 
known distances; solid lines - the upper and lower critical limits M tr 
: (1) M crit for Mi = 1.0, Y = 1.0; (2) M^, for Mi = 0.6, Y = 1.0; 
(3) M c + rit for M, = 1.0, Y = 0.96, Z = 0.04; (4) M+ nt for Mi = 1.0, 
Y = 1.0; (5) M+ it for Mi = 0.6, Y = 1.0. (Mi - in sotar units). Dotted 
lines: Evolution models for AM CVns through the WD channel (kindly 
provided by Chris Deloye). 



The critical rate M CT for helium is 12 times higher than for 
hydrogen-dominated discs consequently the existence of cold 
stable hydrogen-deficient discs does not require ridiculously low 
mass-transfer rates as in the case of hydrogen-rich CVs, where 
it has to be lower than 8 x 10 12 g/s. The outer disc radius, Ru, is 
determined from the equation 



Rd 

a 



0.6 



(1+q) 



1/3' 



(15) 



(Warner 2003), where a is the binary separation and q = Mi/M\ 
the mass ratio. Since in AM CVn stars 0.0125 < q < 0.18 the 
dependence of the maximal disc radius on mass ratio is rather 
weak and can be safely neglected. 

The inner radius, R m , is taken to be equal to the radius of 
the central white dwarf, = Rwd, which is determi ned by the 
white dwarf mass Mwd through the M - R relation dNaue nberg 
|1972|) . The dependence on a in equations ( [TBI and (fT4l is clearly 
negligible. The plot of the estimated M tr versus P 01 -b for observed 
AM CVns (see Table 13 compared with the model predictions is 
shown in Figj3] In agreement with the model, the permanently 
high-state systems ES Cet, SDSS J1908+3940 and AM CVn lie 
well above the critical mass transfer limit. Another bright and 
steady system is HP Lib. Although it is very close to the crit- 
ical line of M + .,, there still exists a certain set of M\ - lines 

cnt' cnt 

above which HP Lib lies with reasonable parameters. The mass- 
transfer rates of CR Boo and V803 Cen are very close to the 
upper critical line, which might explain why these two outburst- 
ing systems have been alternately and confusedly classified as 
analogues of ER UMa, Z Cam, and (not consistent with their 
mass-transfer rate) of WZ Sge stars. Below we test the hypoth- 
esis that they are 'superoutbursting" Z Cam-type stars. GP Com 
and V396 Hya have shown no changes in their brightness till 
now, they are considered to be low-state systems. From Eq. (fT4l 
one concludes that their primary masses Mi should lie in the 
range 0.6 - 1 .0 M Q (the mass-transfer rate for GP Com is an up- 
per limit, see Table 0. 
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System 


r orb (mm; 


Mi(NIq) 




iOg M lr Lg/SJ 


tvcl. 


ES Cet* 


10.35 


0.44 - 0.69 


0.062 - 0.26 


1 *7 TT+O.OV 

1 '-^'-0.19 


1,5 


SDSS J1908+3940* 


~ 15.6 b 


7 


9 


17.35 - 17.73 


7 


AM CVn* 


17.15 


0.71 ±0.07 


0.13 +0.01 


17.65+°'[ 2 


2 


HP Lib 


18.38 


0.80 - 0.4 


0.048 - 0.088 


16 98 +017 
lu.^o _ 027 


2 


CRBoo 


24.52 


1.10-0.67 


0.048 - 0.088 


16 7 +0 18 

lu - -0.32 


2 


V803 Cen 


26.6 


1.17 - 0.78 


0.059 - 0.109 




2 


SDSS J0926+3624 


28.3 


0.85 ± 0.04 


0.035 ± 0.003 


15.79+U, 


4, 6 


SDSS J1240-01 


37.36 


? 


? 


14 29+ 032 


3 


GP Com* 


46.57 


0.50-0.68 


0.009-0.012 


< 14.35^;f 


2 


V396 Hya* 


65.1 


? 


? 


13.33 


3 



Table 2. Prope rties of AM CVns wi th kn own orbital pe r iods and distan ces. *- persistent s yste ms; t - upper l i mit, b - sup erhump period (?). 



Referen ces: CD lEspaillat et all d2005h : (2) iRoelofs et al.N2007al) : (3) iBildsten et all d2006h : (4) iDelove et al] d2007h : (5) iCo pperwheat et al] 



d2011bT) : (6)[Copperwh eat et al. 



fcOllal) . (7) lFontaine et al.l 420111) . See also iNelemansf d20 1 Oh . 



5. Modelling outbursts of helium-dominated 
accretion discs 

5. 1 . The front structure 

The structure of heating and cooling fronts is affected by the 
atomic structure of the dominant elements. Helium has two ion- 
ization energy levels: £Heii = 24.6 eV which corresponds to 
temperature Tu e u ~ 28500 K for the ionization of the first elec- 
tron, and Eueiu = 54.4 eV, i.e. 7H e m ~ 63000 K, for the ioniza- 
tion of the second electron. The existence of two electrons in the 
helium atom also opens a wide range of states where one elec- 
tron is exited and the second one is ionized, or both are excited 
to various energy levels. 

For solar composition discs, a detailed analysi s of the cool- 
ing an d heating fronts has been performed by iMenou et al.l 
(1999). They found that the form of the fronts does not de- 
pend on viscosity and front location in the disc. In hydrogen- 
dominated discs the temperature gradients of both heating and 
cooling front are quite smooth (the bottom row of Fig. |4]shows 
surface density (left) and central temperature (right) profiles 
in solar-composition disc), showing no clear breaks. In com- 
parison, fronts in helium discs exhibit additional features (top 
row of Figj4]i. The breaks in the slope of the temperature gra- 
dient appear at roughly the same temperatures for the heat- 
ing and cooling fronts. One can identify these temperatures 
with atomic transitions in helium. The first break appears for 
T = 80000 - 92000 K (69 - 80 eV) and can be assigned to 
double-photoionization and different negative-ion resonances. 
The break at T ~ 50000 K (~ 42 eV) corresponds to the state 
in which the atom is singly ionized and the second electron is 
in excited state n = 2. In the region of temperatures around 
T ~ 28000 K (~ 24 eV) the neutral helium becomes singly ion- 
ized. Below this is a range of temperatures at which one electron 
remains in the ground state while second one is in the excited 
state. Other points marked on T and X profiles are hard to clas- 
sify unambiguously. 

Yet more changes with the addition of metals to helium. We 
consider here discs with Y = 0.98, Z = 0.02 and Y = 0.96, 
Z = 0.04. At the temperatures available in these discs, elements 
such as oxygen or carbon become highly ionized. Despite low 
number fractions, their contribution to the population of free 
electrons, and in turn their impact on the opacities in the disc, 




0.35 0.4 0.45 0.5 0.55 0.35 0.4 0.45 0.5 0.55 

R[10 1D cm] R[10 ,0 om] 




0.3 0.35 0.4 0.45 0.5 0.3 0.35 0.4 0.45 0.5 

R[10 ,0 cm] R[10 10 cm] 



Figure 4. Cooling front structure: top, left: Y=l, £ profile; top, right: 
Y=l, T c profile; bottom, left: solar, E profile; bottom, right: solar, T c 
profile. Parameters used for calculation are a c = 0.05, 07, = 0.2, Mi = 
1.0, M a = 1.0 x 10' 5 g/s, (R D ) = 1.2 x 10'° cm. 



is significant. For instance, the fourth ionization energy of O is 
Eov = 77.7 eV (T ~ 82000 K), the fourth ionization energy of 
C is E C v = 64.7 eV (T ~ 68000 K) and that of N (thought to 
be enriched in some AM CVn stars) is £nv = 77.47 eV (T ~ 
81760K). Because the main opacity sources in the disc are free- 
free and bound-free transitions (electron scattering, however, is 
negligible), the additional electrons change the efficiency of the 
cooling and heating mechanisms. This influences the maximum 
and minimum temperatures in the disc (see Table [TJ. N enrich- 
ment does not contribute more than other metals to the opacity 
growth. For solar-composition discs we do not see the features 
from atomic transitions of helium and metals, because fronts 
there start to propagate at lower temperature and much lower 
densities (i.e. hydrogen ionization temperature ~ 3.5 x 10 4 K 
and E ~ 30 g cm -2 ) at which these transition do not take place 
or are negligible because of the small contribution of metals to 
the overall disc composition (such as single ionization of carbon 
at T ~ 12000 K ). 
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Figure 5. Light curves calculated within the standard DIM with the 
same a in the hot and cold state for the discs with two different chemi- 
cal compositions. Top: Y=1.0. Bottom: Y=0.98, Z=0.02. Parameters in 
both cases : a h = a c = 0.1, M„ = 10 16 g/s, Mi = 1.0 and the mean 
outer radius is <Rd> = 1-2 x 10 10 cm. 



The influence of the chemical composition on the details of 
the front structure does not manifest itself in the observational 
outburst cycle but helps understanding the DIM physics. 

5.2. Playing with a' s 

The critical surface densities depend on a to the power ~ 0.85 so 
that an or's ratio of 4, say, corresponds to Z's ratio of ~ 3.2. From 
Table Q] we can see that this is just the ratio of critical densities 
for a pure helium disc with a/, = a c . Indeed, for such a disc 
the DIM produces outbursts with amplitudes up to ~ 2.5 mag 
with a kept constant (Fig. [5). However, this cannot correspond to 
real AM CVn outbursts since, as mentioned above, observations 
show that their discs also contain certain amount of metals. We 
tested our model for discs with metallicity Z = 0.0003, 0.02, 
and 0.04. For Z > 2% one cannot obtain amplitudes larger 
than 1 mag without changing a. Adding just 2% of metals low- 
ers the amplitude to less than 1 mag (Fig. 0. Nevertheless, as 
we will below, in AM CVn the required a jump can be lower 
than in dwarf novae (see e.g. Fig. |7). One can notice also that 
changing the chemical composition modifies not only the am- 
plitude but also the shape of the outburst light-curve. This is 
to be expected because a diminished E~ .. ratio favours the 

1 cnt' cnt 1| 1 

appearance of so-called "reflare s" (see e.g. IMenou et al.l 12000: 
iDubus et al.ll200UlLasotall2001h . 

The accretion rate at outburst maximum is M pea k ~ 
{ati/acf^ 1 (see e.g. IDubus et alJl200ll) . hence the necessity of 



increasing the ratio of the viscosity parameters. There is not 
much to guide us in performing the change of a. Numerical 
simulations of the magneto-rotational instability (MRI) which 
is believed to drive turbulence in Keplerian accretion discs (see 
Balbusl l2003L 12005, and references therein), provide no guide 
in this case. Indeed, the value of a they produce is an order of 
magn itude lower than required by observations (e.g. Kin g et all 
1200% . Transients appearing in some MRI simulations produce a 
viscosity coefficient of the correct order of magnitude but they 
are too short-lived to be the solution of this fundamental problem 
(ISorathia et al.ll20fll and Reynolds, private communication). 

The choice of the formula matching the two as is not without 
consequences becaus e it can affect the prop erties of the outburst 
light-curves (see e.g. lCannizzo et al. 2010). In our code we use 
the following interpolation formula: 



log a(T c ) = log a c + [log a h - loga c ] x 



1 + 



(16) 



The motivation behind this expression is the wish to keep the 
values of I^(a) and 2 + (a) in the effective S curve equal to the 
"original" and (H98). Much depends, however, on the 

*^ cm cut A I | I i 

choice of T . H98 chose T « 7/ c (£ crit ) while lDubus etalJ (12001) 
used Tq = 0.5[r c (E~. t ) + T c (T,t. )], which is more convenient for 
the irradiated disc they were considering (irradiation changes the 
critical values on the hot branch). In the present paper we are 
using the H98 prescription. 

5.3. Additional sources of disc heating 

Additional sources of disc heating are known to be important in 
H-rich CVs, and we need to test their influence on the modelled 
light-curves for helium discs. 



5.3.1 . Outer disc heating 

As the mass is transferred from the secondary, its stream hits the 
outer ring of the accretion disc, forming the so-called hot spot. 
In the impact region the dis c annulus of width AR^s is heated at 
a rate Qi with efficiency rjj (Buat-Menar d et al . 2002): 



Qi(R) = n 



GM\M tr 



1 



2R D 2ttR d AR, 



-exp - 



Rn-R 



ARh 



(17) 



where M\ is the mass of the primary white dwarf, M tr is the mass 
transfer rate from the secondary and Rd is the outer radius of the 
disc. 

When the system is in the high state most of its luminos- 
ity comes from the hot accretion disc, especially from the inner 
parts through dissipation in the boundary layer, but during the 
low state as much as ha lf of the visible light can come from 
the hot spot dSmakll2bl0h . The magnitude of this contribution 
depends on the inclination of the disc: lower inclination means 
higher luminosity from the impact region. 

When the term Qi(R) is taken into account in energy con- 
servation equation, the heating from the hot spot reduces the 
critical values of S, which facilitates triggering the outside- 
in outbursts for mass transfer rates lower than in absence of 
heatings. It also decreases the lower limit of the M tl range for 
which the disc should be unstable, and lifts the quiescence level. 
In several cases including this term in the equations improves 
the description of dwarf nova outbursts consider ably (see e.g . 
iBuat-Menard et alj 1200 llfbh . On the other hand, ISmakl (120021) 
expressed doubts concerning the relevance of hot-spot heating, 



7 



Iwona Kotko et al.: Models of AM CVn star outbursts 



noticing that most of the impact energy is radiated away at the 
hot spot location. 

Another presumed source of the outer disc heating i s the ac- 
tion of tidal stres ses (see e.g. iBuat-Menard et al . 2001). In this 
case ISmakl (12001 argues that effect is important only during out- 
burst and limited strictly to the outer disc edge. 

The two outer disc heating contributions are of the same or- 
der of magnitude. In the following, when testing the DIM's prop- 
erties, we will use only hot-spot heating (EgfTTTi because it rep- 
resents the presumed effect on the outer disc structure. 

5.3.2. Irradiation by the hot white dwarf 

As shown bv lHameurvet all (119991) . in dwarf novae irradiation 
by the hot white dwarf can seriously affect outburst structure and 
evolution. In addition to stabilizing the innermost disc regions by 
maintaining their temperature above the ionization temperature, 
this irradiation has a destabilizing effect on the adjacent parts 
of the flow. There, irradiation reduces the vertical disc tempera- 
ture gradient which decreases the critical surface density E~ ri . 
In hydrogen-dominated discs T~ ~ 8000K while the accret- 
ing wh ite dwarf temperature reaches values from 15 000 to 50 
000 K dSion & Godonl l2007). Effective temperatures of accret- 
ing w hite dwarfs in AM C Vns are expected to be i n the same 
range (Bildsten et al. 2006), and apparently they are dSion et al.l 
2011). However, the critical midplane temperatures of helium- 
dominated discs are in the range T~ ~ 18000 - 25000K which, 
in general, reduces white dwarf irradiation to an unimportant ef- 
fect that will be neglected below. 

5.4. Magnetic field of the primary 

Observations of SDSS J080449.49+ 161624.8, whose optical 
spectra have properties similar to those of the hydrogen-rich in- 
termediate polars, suggest that th e primary's magneti c field also 
plays a role in AM CVn stars dRoelofs et al .1 12009). If this is 
the case, the white dwarf magnetic moment will have a non- 
negligible influence on the inner part of the disc. The magnetic 
pressure increases steeply with decreasing radius (P mag ~ R b ) 
and in quiescence it might exceed the gas and ram pressures of 
the infalling matter up to radius Rm, disrupting the matter flow. 
During the outburst the mass accretion rate sharply rises and the 
ram pressure of matter dominates over the magnetic pressure so 
thatthe inner edge of the disc approaches the surface of the white 
dwarf. This me chanism is summa rized in the formula for the in- 
ner disc radius dFrank et"aT1l2002l) : 

R in = R M = 9.8 x lO & M^ ,7 M~ in fi^cm, (18) 

where ^30 is the magnetic moment in units of 10 30 Gcm 3 , M\ 
is the mass of the primary in solar masses, and M\s is the mass 
accretion rate in units of 10 15 g s . 

Including formula (fT8l in the model prolongs the quiescence 
time and outburst duration. The explanation is the same for small 
Mi : when R m is larger, more mass has to be accumulated to cross 
^ZjXR), therefore, the quiescence time lengthens. But because 
more mass is now in the disc, its accretion and the decay from the 
outburst will take longer. The magnetic field does not influence 
the outburst amplitude. 

6. Normal outbursts of AM CVn stars 

From the theoretical point of view, normal outbursts of accretion 
discs are best defined as outbursts in which only the thermal- 



viscous instability is at play with no significant mass-transfer 
variations an d no decisive role of the tidal interactions (see 
lLasotall2001l for a detailed description of normal outbursts). 

As mentioned before, there exists a clear difficulty with test- 
ing the (helium) DIM on AM CVn outbursts since, in these sys- 
tems, most of th e time it is not c l ear when and if an observed out- 
burst is normal. [Patterson et al.l (119971 12000) proposed using the 
empirical relation between the observed outburst amplitude and 
the recurrence time (the "Kukarkin-Parenago relation"; hereafter 
"KP relation") as a test of "normalcy" of outbursts. To evaluate 
the relevance of this suggestion, we will now briefly discuss this 
empirical amplitude vs recurrence-time rel ation and its connec- 
tion to the DIM (see Kotko & Lasota 2012, for more details.) 

6.1. The Kukarkin-Parenago relation for AM CVn stars 

Kukarkin & Pa renagol dl934l) suggested the existence of a rela- 
tion between the amplitudes of cataclysmic-variabl e star out- 
bursts and their recurrence times. As clarified by Ivan Par adiis 
(Il985l) this correlation may indeed to be considered to exist if the 
sample considered is reduced to normal dwarf nova outbursts. 
Because of the large scatter of the observational data the form of 
the correlation is sample-dependent and uncertainties are rather 
large. For a selected subset of U Gem type dwarf novae with 
well-measured amplitudes and recurrence times (see Table IB. II ) 
the correlation has the forrrfl 

A„ = (1.3 ± 0.6) + (1.6 ±0.3) log T n , (19) 
where T n is measured in days and A n in magnitudes. 




log T n [days] 

Figure 6. Amplitudes and recurrence times of normal dwarf nova out- 
bursts and the Kukarkin-Parenago relation, triangles: U Gem stars; 
square: PTF1 J07 19+485 8 (normal outburst); Cir cles: V803 Cen and 
CR Boo (cycling state rebrightenings). Data from Warner (2003]) and 
iRitter & Ko lb (2003). The dotted lines represent the upper and lower 
uncertainty of the observational KP relation fitted to the systems marked 
on the plot while the two solid lines represent the theoretical KP rela- 
tion for helium discs - the lower line has been calculated for a h = 0.2, 
R D = 0.7 x 10 10 cm, Mi = 1.0 and the upper line for a h = 0.2, 
R D = 1.0 X 10 10 cm, Mi = 1.2. 

The only AM CVn system with unambiguous normal out- 
bursts is PTF 1J0719+4858. Their amplitude is A n ~ 2.5 and 

3 For the sample in lWarnerl (120031) the relation is A n = (0.7 ± 0.43) + 
(1.9 ±0.22) log T n . 
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recurrence time T„ ~ 10 d. Fig. [6] shows the KP relation 
and the corresponding amplitudes and recurrence times of se- 
lected U Gem stars (marked by triangles). The dotted lines rep- 
resent the upper and lower uncertainty of the KP relation (Eq. 
|T9l >. Also plotted are the amplitudes and recurrence times for 
PTF 1J0719+4858 and the cycling states of the AM CVn stars 
V803 Cen and CR Boo. PTF 1J0719+4 858 fits the KP relation 
well. The same is true, as remarked by [Patterson et al.l (l!997[ 
2000), of the parameters of the cycling states of the two other 
AM CVns: they are compatible with this relation. However, does 
thos imply that they can be classified as normal outbursts? This 

depends on what the KP relation is supposed to r epresent. 

Us ing the standard version of the DIM, iKotko & Lasotal 
(120 121) derived a theoretical A n (T n ) relation in the form of 

A n «C, +2.5 log T n , (20) 

with 

Ci « 1.5 - 2.5 log r dec + BC mdX - BC min , (21) 

where fd ec is (as before) the outburst decay time calculated using 
the disc size and the viscous speed (see Eq. Q]) and BC mdx and 
BC m i n are the disc bolometric corrections in the high and low 
states, respectively. The viscosity parameter a/, must be in the 
range 0. 1 - 0.2 to correspond to observations. 

The relation Eqs. (1201 . ( |2T1 i is based on the assumption that 
the whole mass accumulated during the time f quiesc + fese is ac- 
creted onto the white dwarf during outburst (see Kotk o & Lasotal 
1201 2l for details). Normal outbursts are supposed to satisfy this 
assumption therefore eruptions following the theoretical A n (T n ) 
relation can be considered to belong to this category. In Fig. 
[6] we plotted two theoretical A n {T„) relations corresponding to 
helium discs with parameters Rd = 0.7 x 10 ll) cm, M\ = 1.0 
(PTF 1J0719+48583) and R D = 1.0 x 10 ll) cm, M, = 1.2 
(V803 Cen and CR Boo), respectively. In both cases »/, = 0.2. 
With these parameters the amplitudes and recurrence times of 
helium-dominated discs lie on the theoretical A n (T n ) relation. 

6.2. Model light-curves of normal outbursts 



9.5 




12.5 1 1 1 1 1 1 1 

70 75 80 85 90 95 100 

Time [days] 

Figure 7. Normal outburst cycle for Y=0.98, Z=0.02 accretion disc. 
or c = 0.05, ff,, = 0.1, Mi6 = 0.1, Mi = 1.0, (Rd) = 1.1 x 10 10 cm. The 
disc extends down to the white dwarf surface. No irradiation of disc or 
companion is taken into account. 

4 The orbital parameters of this system are not yet known, however. 



In SU UMa stars the outbursts appearing between the su- 
peroutbursts are believed to be of normal type. The case of 
PTF 1J07 19+4858, where a similar sequence is present, suggests 
that the same is probably true for AM CVn stars. 

Fig. [7] shows a synthetic light curve produced with the DIM 
applied to the helium-dominated discs. Although one does not 
expect to observe the light curve analogous to FigJ7]among com- 
pact helium binaries, the basic DIM should describe this part of 
the outburst cycle where normal outbursts are present. 

We briefly recall how the model's free parameters a/, , a c and 
mass transfer rate M tr influence the recurrence times T„ and am- 
plitudes A n of normal outbursts. The primary mass M\ , fixing 
the inner disc radius R\ n , is often difficult to estimate from obser- 
vations and then it can be considered a free model parameter as 
well. The disc radius is constrained by observations, so there is 
not much freedom left in setting Rd in the model, as long as one 
intends to refer to the real systems. 

The time elapsing between the onsets of two consecutive out- 
bursts extends with increasing ratio B of viscosity parameters: 
B = a/,/a c . Since the fiducial value of a/, is 0.1 - 0.2, it is a c 
which decides about T„. However, B also controls the value of 
A n . Accordingly, the increase of B (by lowering a c ) not only re- 
sults in longer recurrence times, but also in higher amplitudes. 
The outbursts will be more frequent in the disc fed with matter 
at high rates or/and in a system with heavier primary. 

High M t r is also responsible for the appearance of wide nor- 
mal outbur st, whi ch should not be confused with superoutbursts 
(see lLasotall200lh . 

Adding metals to helium increases El^,, due to the additional 
non-negligible contribution from highly ionized metals to the 
opacity. Furthermore, the consequence of augmented Z is the 
reduction of the total amount of helium and E"|\ drop. The con- 
clusion is that adding elements heavier than He to the disc matter 
is equivalent to diminishing the £ ratio (decreasing the B ratio). 
This results in shorter recurrence times and lower outburst am- 
plitudes. 

It is worth emphasizing that metallicity affects the value of 
a c , which needs to be "adjusted" to obtain the observed outbursts 
amplitudes and recurrence times. The difference between a/ t and 
a c to be explained by MRI may not be as large in a helium disc 
as usually thought for a solar-composition disc. However, a dif- 
ference is still required, as argued in Sect. 15.21 



7. AM CVn outburst supercycles 

Light curves of AM CVn can be divided into three categories, 
which we refer to as: KL Dra - like, PTF 1J07 19+4858 - like, and 
CR Boo/V803 Cen - like. The first category exhibits superout- 
bursts only and might be the equivalent of WZ Sge-type dwarf 
nova starsj. PTF 1J0719+4858 almost certainly corresponds to 
the SU UMa type, well-known in hydrogen-dominated dwarf 
novae. The most mysterious are CR Boo and V803 Cen, which 
have been variedly described as equivalent to WZ Sge, SU UMa, 
ER UMa, and Z Cam stars. It seems that they belong to all these 
classes partaking freely of such a wide range of dwarf nova be- 
haviour, to quote Pat terson et al.1 (12000). 

As mentioned in Sect. 12.11 in the present context we have 
chosen to opt for the EMT model as the model explaining the 
mechanism of superoutbursts appearance. 



5 KL Dra might not be really of this type since small amplitude nor- 
mal could have been missed, but an AM CVn analog of the WZ Sge- 
type presumably exists. 
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Although there is ample evidence for both the secondary's 
irradiation an d the enhanced (during outburst) mass-transfer rate 
dSma k 2011), the connection between th e two is far from be- 
ing clear. Viallet & Hameurv (2007, 2008) have shown that mat- 
ter heated up by irradiation of the companion's disc-unshaded 
surface has cooled down too much by the time it arrives in 
the vicinity of the L\ point to affect mass transfer towar ds the 
disc. Although the treatment of the cooling is simplified (Smak 
2009d), it is unlikely that the correct treatment would modify the 
result significantly. A solution for this difficulty would be to have 
direct irradiati on of the L\ region allowed by a tilted/warped disc 
(ISmakll2009et) . 

Below, inspired by previous work, we will attempt to mimic 
the effect of irradiation on the mass-transfer rate by suitable 
parametrizations. 

7.1. KL Dra and PTF 1J071 9+4858 - like outbursts 

The general properties of the relatively simple light curves of 
KL Dra and PTF 1J0719+4858 can be reproduced by using the 
simple prescription for accreti on-irradiation mass-tra n sfer rate 
increase originally proposed by Hameu rv et al.l(ll997l) : 

M tt = max(M ,tr, yM acc ), (22) 

where < y < 1. Mo.tr corresponds to the "secular" (non- 
enhanced) mass-transfer rate while M acc is the accretion rate 
onto the primary. 

Dealing with better-sam pled multi-wavelengt h light curves 
such as that of VY Hyi , Schreiber et al.l (120041) used a more 
"refined" prescription using a suitably averaged over time M acc 
but using it here would be exaggerated. 

As for hydrogen-dominated dwarf novae (lHameurv et al.l 

2000) the light curve properties depend on the mass-transfer rate, 
the mass of the primary, and the viscosity parameters, as well as 
on the assumed white dwarf's magnetic moment and the value 
of the parameter y. In AM CVns one has additionally to fix the 
chemical composition, i.e. the metallicity which we will assume 
to be Z = 0.02 (Y = 0.98) unless stated otherwise. We assumed 
(Rd) = 1.0 x 10 10 cm for the average size of the accretion disc. 

7.1.1. Superoutbursts only 

To obtain a light curve similar to t hat of KL Dra, we us ed param- 
eters close to those suggested by Ramsav et alJ (I2010I) for this 
binary:M tr = 2x 10 16 gs _I , M x = 0.6M o . With the choice of the 
viscosity parameters a c = 0.035, 07, = 0.2, choosing y = 0.8 and 
using = 1 (corresponding to a magnetic field B ~ 1.5xl0 3 G) 
we obtained the supercycle shown in Fig. [8] The calculated am- 
plitude A s « 3 mag, the (super)outburst duration T^ur ~ 15 d 
and the recurrence time r recc » 57 d correspond quite well to the 
observed A s « 3.5 mag, T^m- ~ 14 d and r reC c ~ 63 d. 

Including the magnetic field of the primary and the heating 
by the hot spot was essential for obtaining the light curve in Fig. 
[8] These two effects stabilize the disc. The magnetic field trun- 
cates the inner parts of the disc, increasing 2"" ri .(/?i n ). The disc 
then has to accumulate more mass to trigger the inside-out out- 
burst. The hot spot heating, in turn, lowers Ei^CRo) so that more 
mass has to be accreted before the cooling front can form. This 
favours triggering of a superoutburst instead of a normal out- 
burst. Both effects contribute to the suppression of normal out- 
bursts. The resulting light curve compares well with the obser- 
vations (Fig. [HJ. 

In reality, as mentioned above, KL Dra could have normal 
outbursts that have been missed during observational campaigns. 
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Figure 8. Top : KL Dra model: a c = 0.035, a h = 0.2, M tt = 2 x 
10 16 g/s, Mj = 0.6, y = 0.8, <R D > = 1.2 x 10 10 cm, ^30 = 1. Bottom: 
The part of KL Dra light curve showing two subsequent superoutbursts. 
The time between their maxima is 60 days, their amplitude is 3 mag 
and their duration is ~ 13 days. Data provided by G. Ramsay. 



If true, this would mean that one should drop from the disc evo- 
lution equations e.g. the hot-spot heating term to obtain normal 
outbursts. This is what we have done in the next subsection when 
calculating models of the AM CVn equivalents of SU UMa stars. 

7.1.2. Helium "SU UMa stars" 

Except for its orbital period (»27.7 min), not much is known 
about the parameters of the first observed helium SU UMa-type 
system PTF 1J07 19+4858 . We assumed that this type of system 
transfers mass close to its secular mean taking for example M tr = 
6 x 10 16 gs-'. Assuming Mi = 1.0 M , a c = 0.02, a h = 0.1, 
y = 0.6 and a disc extending down to the white dwarf's surface, 
we obtain the light curve presented in Fig. [9] 

Between superoutbursts one notices a slight gradual increase 
of the minimum and maximum luminosity of the four consec- 
utive normal outbursts. This happens because during each nor- 
mal outburst M tI is enhanced due to the y-prescription - the disc 
gains more mass than it looses due to accretion. In consequence, 
during the sequence of the normal outbursts, the mass is gradu- 
ally accumulated in the disc and S rises everywhere in the disc. 
This effect is one of the typical features of the EMT model (it is 
also present for slightly different reasons in the TTI version, see 
Schreiber et al .1 120041 iTsugawa & Osakilll997h . While the gen- 
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Figure 9. SU UMa type model for a helium-dominated disc (an "imi- 
tation" of PTF 1J0719+4858 ): a c = 0.02, a h = 0.1, M„ = 6 x 10 16 g/s, 
Mi = 1.0, y = 0.6, (R D > = 1.0 x 10'° cm, fi 30 = 0. 



eral luminosity rise in quiescence is not observed in real sy stems 
and i s one of the well-known weaknesses of the DIM (Smak 
2000), the increase of the normal outburst peak luminosity might 
have been observed in some systems. 

The supe routburst precursor is also typical of the EMT and 
TTI models ( Schre iber et al.l 120041) . Owing to mass accumula- 
tion, the disc arrives to a state where after the rise of a normal 
outburst the cooling front is no longer able to propagate - the 
disc becomes stuck in the hot state and a superoutburst begins. 
This last normal outburst leading to a superoutburst appears as 
its precursor in the light curve. 

When one compares the model light curve to that of 
PTF 1J07 19+4858 one notes that in our model the amplitude of 
normal outbursts is larger by 0.5 mag and their duration is four 
times longer. The normal outburst duration is extended by the 
M tl - enhancement and the high amplitude is the result of the as- 
sumed a ratio - a lower ratio would result in a lower amplitude. 
However, lowering this ratio would also lower the superoutburst 
amplitude and shorten the recurrence times. Tuning the param- 
eters to obtain a better agreement does not make much sense in 
view of the arbitrariness of the mass-transfer prescription and of 
the uncertainties of the DI M itself (for a discussion of the weak- 
nesses of the DIM see e.g. Lasota 2001). 

7.2. Dips, "cycling", and standstills 

Most (maybe all) non-stationary AM CVn stars exhibit at least 
one "dip" during the decay from the superoutburst maximum. In 
some systems, such as CR Boo and V803 Cen, the system be- 
comes stuck in a state in which series of ~ 1 .0 mag amplitude, 
~ 0.8 - 1 .0 d recurrence-time outbursts are present. There is also 
some uncertainty about the nature o f these intermissions i n the 
superoutburst cycle: on the one hand. lPatterson et all (1 19971) con- 
sidered them to be a "cycling state" o f normal dwarf nova out- 
bursts, on the other, Kat o et al.1 d2001al) interpreted the same part 
of the (CR Boo) light curve as a standstill (analogous to that of 
Z Cam dwarf novae), arguing that no superoutbursts are present 
during the high-luminosity state maintained by the system. 

To answer which of these two interpretation is correct more 
observations are needed and therefore we are not trying here to 
decide this question. Let us just note that the mass-transfer rates 
attributed to both CR Boo and V803 Cen (see Table |2]i are very 
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Figure 10. Light-curve of a superoutbursting helium "Z Cam" star. 
The parameters are: a c = 0.04 , a h = 0.2 , M\ = 1.0, Mo.tr = LI x 
10 17 g/s, AM tr = 0.4M 0>tr , T m = 10 d, r low = 70 d, y = 0.7, <R C > = 
l.Ox 10'° cm. The dotted lines correspond to the hot stability limit. Top: 
the model lightcurve. Bottom: M v modulations. 

close to the critical ones and a phenomenon analogous to Z Cam 
standstills could indeed be expected for this systems. The main 
difference is that Z Cam stars ha ve no superoutbursts. We w ill 
therefore adapt the Z Cam model dBuat-Menard et aD 2001b) by 
combining its mass-transfer modulations with the EMT. 

Dips might or might not be related to the cycling states. They 
are very similar to those observed during the decay from super- 
outbursts in WZ Sge - type stars. If common to both systems, the 
dip origins are not connected to the value of the mass-transfer 
rate but would instead result from the very compact size of their 
orbit. We tested a simple hypothesis based on this assumption. 

7.2.1 . "Z Cam - type" modulations and (y) superoutbursts 

Adapting the Z Cam models of iBuat- Menard e tal] (l2001bl) to 
AM CVn stars, we assumed the mass-transfer to be modulated 
as AM /(id) - 40% around an average rate close to the crit- 
ical (Mo.tr = 1.1 x 10 17 gs~') with the modulations occurring 
at t = 10 + k ■ 70 days. To take into account the presence of 
irradiation-induced super-outbursts, we combined this modula- 
tion with the mass-transfer modulations given by Eq. (l22t with 
the actual value of the mass-transfer rate instead of Mo. tr . 

Choosing as typical parameters a c = 0.04, a/, = 0.2, Mi = 
1.0 M , (Rd> = l.Ox 10 10 cm and y = 0.7 we obtained the 
light-curve shown in Fig. [10] This light curve shows superout- 
bursts and standstills, as expected, but features resulting from 
the superposition of both types of modulations are also present 
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in the form of "dips". The hypercycle (the cycle between two 
consecutive superoutbursts followed by a standstill) starts with 
a major M tr enhancement due to y (about day 70 on FigJTOb. 
During the decay from the superoutburst maximum, the fall of 
M tt is quenched by arising modulation through Mo tl +AM tl . The 
mass-transfer rate then rises again (around day 80) and the de- 
cay from maximum is reversed by a resulting outside-in heating 
front that catches up with the propagating cooling front before 
the latter reaches the inner disc edge - hence a dip-like feature. 
During the following superoutburst a slight difference in phase 
of the two mass-transfer modulations allows the cooling front to 
propagate almost to the disc's inner end. A similar mechanism 
produces the dip in narrow outburst preceding the superoutburst. 
This suggest that dips and cycling-state features might result 
from mass-transfer modulations triggering cooling/heating front 
"catching" and not from reflections. On the other hand, the in- 
dentation at the end of the sta ndstill is the result o f cooling front 
reflection (as described e.g. in Dub us et alj|1999l) . 

In any case, the presumed Z Cam-effect mechanism applied 
to outbursting AM CVns produced the required type of light- 
curve combining superoutbursts with standstills (and normal 
outbursts). 

7.2.2. Sinusoidal plus y modulation 

The prescription for mass-transfer modulations used in the pre- 
vious section has been chosen not for its realism but because 
it has been used in modell ing the Z Cam star outburst cycle by 
iBuat-Menard et al.1 d2001bl) . Since in AM CVn light curves other 
peculiar features are present in addition to standstills, we tried 
other simple forms of mass-transfer rate modulations. The sim- 
plest is to modify Eq. (1221 as 

M tI = max(M , tl -( 1 + A sin(C + nt/r), (yM acc ))), (23) 

where A, C (dimensionless) and t (time) are adjustable con- 
stants; t is the time coordinate. 

With Mi = l.OM , M tr = 1.1 x 10 17 g/s and (R D ) = 
1.0 x 10 10 cm, a c = 0.04, a h = 0.2, choosing A = 0.5, C = 2, 
r = 1 d and y = 0.8, one obtains the light-curve shown in Fig. 
[TT1 One remarks superoutbursts with amplitude A s ~ 3 mag, du- 
ration r s ~ 15 d and recurrence time of r reccs ~ 27 d. These are 
"standard" EMT, y enhancement triggered outbursts. 

In addition, short outbursts with a repetition time of ~ 1 .8 d 
and an amplitude ~ 1.2 — 2.4 mag appear during the decay 
from superoutburst and persist until the following one. One is 
of course tempted to identify them with the cycling states of 
CR Boo or V803 Cen. The origin of these short outbursts is quite 
simple to understand. At the end of the decay of M tr from the 
superoutburst value, the sinusoidal enhancement of M tl domi- 
nates over the y-induced M a fluctuations. This short enhance- 
ment of M tr results in a low-amplitude, short outburst (starting 
at day ~ 73 in FigJTTb because the cooling front that switches off 
the superoutburst is quenched by a freshly triggered outside-in 
heating front. This happens because the post cooling-front S is 
still close to the critical value due to M u enhancement through 
the sinusoidal modulation. However, since the disc's mass is de- 
creasing, after one day a new cooling front forms immediately 
and this time it succeeds in propagating till the inner disc's rim 
but, because of mass accumulating, a new outburst starts im- 
mediately. All heating fronts are triggered by mass-transfer en- 
hancement and therefore are of the outside-in type. Their maxi- 
mum brightness is slightly increasing because of mass accumu- 
lation. Their amplitude depends on the combined strengths of 



the sinusoidal plus y enhancement. If the enhancement brings 
M tr above Mj, the resulting short outburst will have a small am- 
plitude as the outside-in heating front forms easily and promptly 
catches up with the cooling front that quenches it. But when the 
disc's mass is still quite low and the y + sine enhancement is not 
powerful enough to bring M tI above M+, the cooling front will 
propagate unhindered till the low-luminosity level (see e.g. day 
81). It seems therefore that suitable modulations of the mass- 
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Figure 11. Light-curve corresponding to sinusoidal + y modulations 
of the mass-transfer rate. The parameters are a h = 0.2, a c = 0.04, 
M ,tr = 1-1 x 10 17 g/s, Mi = 1.0, <R > = 1.0 x 10 10 cm, y = 0.8, 
A = 0.5, C = 2, t = Id. The dotted lines correspond to the hot stability 
limit. Top: the model light-curve. Bottom: M tl modulations. 



transfer rate can produce various observed features of outburst- 
ing AM CV light curves even with the simplest assumptions 
about their shapes. However, the mechanisms that are able to 
produce mass-transfer modulations are still to be firmly iden- 
tified. There is no doubt that mass-transfer rates in close bina- 
ries are highly variable, in some cases with huge amplitudes. 
Although irradiation of the secondary is the usual suspect, in 
many cas es the observe d variability (e.g. of AM Her or VY Scl 
stars, see Warner 2003) is not caused by irradiation. A periodic 
modulation of mas s transfer ra te could be also achieved by a 
warped/tilted disc (Smak1 l20Tll) . The sinusoidal variation of Eq. 
d23l could be considered an attempt to represent the effect of the 
variable irradiation if the modulation time r = Id could find a 
convincing interpretation in this context. 
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8. Discussion and conclusion 

We have presented a systematic study of the dwarf nova disc 
instability model applied to outbursts of AM CVn stars. What 
distinguishes these hydrogen-free compact binaries from their 
cataclysmic variable cousins is (obviously) the chemical com- 
position and the shorter orbital periods. Both can affect the out- 
burst properties of AM CVn stars. However, because of the still 
unknown mechanism driving the in this context necessary vari- 
ations of viscosity (or the angular-momentum and disc-heating 
mechanisms), it is difficult to assess the importance of helium 
richness in the outburst processes. They are potentially impor- 
tant, especially for very low metallicities, but their real impact 
must await a better understanding of the basic mechanisms driv- 
ing accretion in discs in general. 

The main conclusion of our study is that outbursts of 
AM CVn stars can be explained by the suitably adapted dwarf- 
nova disc instability model. As in the case of the application 
of this model to hydrogen-dominated cataclysmic variables, one 
has to resort to additional effects to account for the observed su- 
peroutbursts, dips, cycling states, and standstills. The basic fea- 
ture that must be included in the model is the variable mass- 
transfer rate. In particular, we have shown that the enhanced 
mass-transfer rate, due presumably to variable irradiation of 
the secondary, must not only be taken into account but is a 
factor that determines the shape of most AM CVn star out- 
burst light-curves. The cause of variable secondary's irradiation 
has yet to be understood; the best candidate is the precession 
of a tilted/warped disc. We postulate the existence of super- 
outbursting Z Cam type A M CVn stars that have been tentatively 
identified in observations dKato et alJl2001h . 

The very short periods and low mass-ratios of AM C Vns im- 
ply the absence of a pure normal-outburst cycle - the test bench 
of the model. The compactness of AM CVn stars is clearly im- 
portant because it accentuates the interactions between the bi- 
nary components. Their mass ratios 0.0125 < q < 0.18 sug- 
gest the importance of resonances and tidal effects in general. 
Irradiation is important as well. Unfortunately, as in the case of 
cataclysmic variables, the implications of these effects on the 
outbursts is uncertain and subject to controversy. At best, they 
can be parameterized only in a fairly rough way. The only effect 
that, presumably important in CVs, can be neglected in AM CVn 
stars is the irradiation of the inner disc regions by the hot accret- 
ing white dwarfs. 

Another drawback that might soon be overcome is the usu- 
ally unsatisfactory quality of observations of the AM CV out- 
bursts. When better and richer sets of data will become available 
the comparison of the DIM with observations will become a pre- 
cious source of knowledge about accretion disc physics. 

Acknowledgements . We are grateful to Chris Deloye for providing us with evo- 
lutionary tracks for AM CVn stars and to Joe Patterson for the CR Boo and 
V803 Cen light curves and for helpful comments, to Gavin Ramsay for the 
light curve of KL Dra and to Eran Ofek and David Levitan for very help- 
ful discussions about the AM CVn stars observed by the Palomar Transient 
Factory. We thank the anonymous referee for his report that helped improving 
the article. This work has been supported in part by the Polish MNiSW grants 
PSP/K/PBP/000392 and N N203 380336, the French Space Agency CNES, and 
the European Commission via contract ERC-StG-20091 1. 

References 

Balbus, S. A. 2003, ARA&A, 41, 555 

Balbus, S. A. 2005, in The Astrophysics of Cataclysmic Variables and Related 
Objects , 330, 185 

Bildsten, L., Townsley, D. M., Deloye, C. J., & Nelemans, G. 2006, ApJ, 640, 
466 



Buat-Menard, V, & Hameury, J.-M. 2002, A&A, 386, 891 
Buat-Menard, V, Hameury, J.-M., & Lasota, J.-P. 2001, A&A, 366, 612 
Buat-Menard, V, Hameury, J.-M., & Lasota, J.-P. 2001, A&A, 369, 925 
Cannizzo, J. K. 1984, Nature, 311, 443 
Cannizzo, J. K., & Wheeler, J. C. 1984, ApJS, 55, 367 

Cannizzo, J. K., Still, M. D., Howell, S. B., Wood, M. A., & Smale, A. P. 2010, 

ApJ, 725, 1393 
Copperwheat, C. M., et al. 2011, MNRAS, 410, 1113 
Copperwheat, C. M., et al. 2011, MNRAS, 429 

Deloye, C. J., Taam, R. E., Winisdoerffer, C, & Chabrier, G. 2007, MNRAS, 
381, 525 

Dubus, G, Lasota, J.-P, Hameury, J.-M., & Charles, P. 1999, MNRAS, 303, 139 
Dubus, G, Hameury, J.-M., & Lasota, J.-P. 2001, A&A, 373, 251 
Espaillat, C, Patterson, J., Warner, B., & Woudt, P. 2005, PASP, 117, 189 
El-Khoury, W., & Li, J. 1999, Annapolis Workshop on Magnetic Cataclysmic 

Variables, 157, 393 
Fontaine, G, et al. 2011, ApJ, 726, 92 

Frank, J., King, A., & Raine, D. J. 2002, Accretion Power in Astrophysics, 
by Juhan Frank and Andrew King and Derek Raine, pp. 398. ISBN 
0521620538. Cambridge, UK: Cambridge University Press, February 2002. 

Hameury, J.-M., Lasota, J.-P., & Hure, J.-M. 1997, MNRAS, 287, 937 

Hameury, J.-M., & Lasota, J.-P. 2005, A&A, 443, 283 

Hameury, J.-M., Lasota, J.-P., & Warner, B. 2000, A&A, 353, 244 

Hameury, J.-M., Menou, K, Dubus, G, Lasota, J.-P., & Hure, J.-M. 1998, 
MNRAS, 298, 1048 

Hameury, J.-M., Lasota, J.-P., & Dubus, G. 1999, MNRAS, 303, 39 

Ichikawa, S., & Osaki, Y. 1992, PASJ, 44, 15 

Idan, I., Lasota, J. P., Hameury, J.-M., & Shaviv, G. 1999, Phys. Rep., 311, 213 
Idan, I., Lasota. J.-P., Hameury, J.-M., & Shaviv, G. 2010, A&A, 519, Al 17 
Kato, T, Nogami, D., Baba, H., Hanson, G, & Poyner, G. 2000, MNRAS, 315, 
140 

Kato, T, Stubbings, R., Monard, B., Pearce, A., & Nelson, P. 2001, Information 

Bulletin on Variable Stars, 5091, 1 
Kato, T, Reszelski, M., Poyner, G, et al. 2001, Information Bulletin on Variable 

Stars, 5120, 1 

Kato, T., Stubbings, R., Monard, B., et al. 2004, PASJ, 56, 89 
King, A. R., Pringle, J. E., & Livio, M. 2007, MNRAS, 376, 1740 
Kotko, I., Lasota, J.-P. 2012, A&A, submitted 

Kotko, I., Lasota, J. -P., & Dubus, G. 2010, Astronomische Nachrichten, 331, 231 
Kukarkin, B. W. & Parenago, P. P., 1934 Var. Star. Bull. 4, 44 
Lasota, J.-P. 2001, New A Rev., 45, 449 

Lasota, J. P., Hameury, J. M., & Hure, J. M. 1995, A&A, 302, L29 

Lasota, J.-P., Dubus, G, & Kruk, K. 2008, A&A, 486, 523 

Levitan, D., Fulton, B. J., Groot, P. J., et al. 201 1, ApJ, 739, 68 

Levitan, D., Kulkarni, S., Prince, T., Nissanke, S., & Vallisneri, M. 2011, APS 

Meeting Abstracts, 1 1008 
Mason, K. O., Cordova, F. A., Watson, M. G. & King, A. R. 1988, MNRAS. 

232, 779 

Menou, K, Hameury, J.-M., & Stehle, R. 1999, MNRAS, 305, 79 
Menou, K., Hameury, J.-M., Lasota, J.-P, & Narayan, R. 2000, MNRAS, 314, 
498 

Nauenberg, M. 1972, ApJ, 175, 417 

Nelemans, G. 2010 AM C Vn stars, GN Wiki 

Nelemans, G. 2005, The Astrophysics of Cataclysmic Variables and Related 
Objects , 330, 27 

Nelemans, G, Yungelson, L. R„ van der Sluys, M. V, & Tout, C. A. 2010, 

MNRAS, 401, 1347 
Osaki, Y. 1985, A&A, 144, 369 
Osaki, Y. 1989, PASJ, 41, 1005 
Osaki, Y. 1996, PASP, 108, 39 
Patterson, J., et al. 1997, PASP, 109, 1100 

Patterson, J., Walker, S., Kemp, J., O'Donoghue, D., Bos, M., & Stubbings, R. 

2000, PASP, 1 12, 625 
Payne-Gaposchkin, C. 1957, The Galactic Novae, North-Holland Publishing 

Company and Interscience Publishers, Chapter 8 
Payne-Gaposchkin, C. 1977, Novae and Related Stars, M. Friedjung ed., 

Dordrecht, D. Reidel Publishing Co. Astrophysics and Space Science Library, 

65,3 

Ramsay, G, Groot, P. J., Marsh, T., Nelemans, G, Steeghs, D.. & Hakala, P. 

2006, A&A, 457, 623 
Ramsay, G, et al. 2010, MNRAS, 407, 1819 

Ramsay, G, Barclay, T., Steeghs, D., et al. 2012, MNRAS, 419, 2836 
Ritter, H., & Kolb, U. 2003, A&A, 404, 301 

Roelofs, G. H. A., Groot, P. J., Benedict, G. E, McAithur, B. E., Steeghs, D., 

Morales-Rueda, L., Marsh, T. R., & Nelemans, G. 2007, ApJ, 666, 1174 
Roelofs, G. H. A., Nelemans, G, & Groot, P. J. 2007, MNRAS, 382, 685 
Roelofs, G. H. A., Groot, P. J., Steeghs, D., et al. 2009, MNRAS, 394, 367 
Schreiber, M. R. 2007, A&A, 466, 1025 



13 



Iwona Kotko et al.: Models of AM CVn star outbursts 



Schreiber, M. R., Hameury, J.-M., & Lasota, J.-P. 2004, A&A, 427, 621 
Shakura, N. I., & Sunyaev, R. A. 1973, A&A, 24, 337 

Shears, J., Brady, S., Koff, R., Goff, W., & Boyd, D. 2011, JBAA, in press, 
arXiv: 1104.0107 

Sion, E. M., & Godon, P. 2007, 15th European Workshop on White Dwarfs, 372, 
563 

Sion, E. M., Linnell, A. P., Godon, P., & Ballouz, R.-L. 201 1, ApJ, 741, 63 

Skillman, D. R., Patterson, J., Kemp, J., et al. 1999, PASP, 111, 1281 

Smak, J. 1983, Acta Astron., 33, 333 

Smak, J. 1984, Acta Astron., 34, 161 

Smak, J. 1984, PASP, 96, 5 

Smak, J. 1989, Acta Astron., 39, 201 

Smak, J. 1999, Acta Astron., 49, 391 

Smak, J. 2000, New A Rev., 44, 171 

Smak, J. 2002, Acta Astron., 52, 263 

Smak, J. 2006, Acta Astron., 56, 365 

Smak, J. 2008, Acta Astron., 58, 55 

Smak, J. 2009, Acta Astron., 59, 89 

Smak, J. 2009, Acta Astron., 59, 103 

Smak, J. 2009, Acta Astron., 59, 109 

Smak, J. 2009, Acta Astron., 59, 121 

Smak, J. 2009, Acta Astron., 59, 419 

Smak, J. 2011, Acta Astron., 61, 59 

Smak, J., & Waagen, E. O. 2004, Acta Astron., 54, 433 

Smak, J. 2010, Acta Astron., 60, 83 

Sorathia, K. A., Reynolds, C. S., Stone, J. M., & Beckwith, K. 2011, 

arXiv:l 106.4019 
Strohmayer, T. E. 2004, ApJ, 614, 358 
Tsugawa, M., & Osaki, Y. 1997, PASJ, 49, 75 
van Paradijs, J. 1985, A&A, 144, 199 
van Paradijs, J. 1983, A&A, 125, L16 
Viallet, M., & Hameury, J.-M. 2007, A&A, 475, 597 
Viallet, M., & Hameury, J.-M. 2008, A&A, 489, 699 
Wade, R. A., Eracleous, M., & Flohic, H. M. L. G. 2007, AJ, 134, 1740 
Warner, B. 1987, MNRAS, 227, 23 

Warner, B. 2003, Cataclysmic Variable Stars, (Cambridge Astrophysics Series, 
28) 

Wood, M. A., Casey, M. J., Garnavich, P. M., & Haag, B. 2002, MNRAS, 334, 
87 



Appendix A: Critical parameters for various 
chemical compositions of helium dominated 
discs. 

I this appendix we summarize the critical parameter formulae obtained for 
various chemical compositions of helium-dominated accretion discs. These 
formulae have been obtained through fits to numerically calculated S-curves. 
One should note that whereas the fit of the critical T c ff is accurate within few 
percent, the accuracy of the fit to the critical central temperature is no better 
than 20%. 
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Appendix B: Amplitudes and recurrence times for a 
selection of U Gem type dwarf novae and 
outbursting AM CVn stars. 





Porb [hr] 


A„ [mag] 


T n [d] 


IP Peg 
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CW Mon 
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1.2 


1 


V803 Cen 


0.44 


1.3 


0.95 


PTF1J0719 


0.45 
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Table B.l. The list of U Gem- type and AM CVn systems (three last 
rows) plotted on Fig. [6] The columns are: P or b is the orbital period, A„ is 
the normal outb urst amp l itude and T n is the norma l out burst reccurence 
time. Data fromjWarner (2003) and Ri tter & KolfJ d2003l) . 
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